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Abstract 

A thought experiment is proposed to demonstrate the existence of a gravitational, vector Aharonov- 
Bohm effect. We begin the analysis starting from four Maxwell-like equations for weak gravitational 
fields interacting with slowly moving matter. A connection is made between the gravitational, vector 
Aharonov-Bohm effect and the principle of local gauge invariance for nonrelativistic quantum matter 
interacting with weak gravitational fields. The compensating vector fields that are necessitated by this 
local gauge principle are shown to be incorporated by the DeWitt minimal coupling rule. The nonrela- 
tivistic Hamiltonian for weak, time-independent fields interacting with quantum matter is then extended 
to time-dependent fields, and applied to problem of the interaction of radiation with macroscopically 
coherent quantum systems, including the problem of gravitational radiation interacting with supercon- 
ductors. But first we examine the interaction of EM radiation with superconductors in a parametric 
oscillator consisting of a superconducting wire placed at the center of a high Q superconducting cavity 
driven by pump microwaves. Some room-temperature data will be presented demonstrating the splitting 
of a single microwave cavity resonance into a spectral doublet due to the insertion of a central wire. 
This would represent an unseparated kind of parametric oscillator, in which the signal and idler waves 
would occupy the same volume of space. We then propose a separated parametric oscillator experiment, 
in which the signal and idler waves are generated in two disjoint regions of space, which are separated 
from each other by means of an impermeable superconducting membrane. We find that the threshold for 
parametric oscillation for EM microwave generation is much lower for the separated configuration than 
the unseparated one, which then leads to an observable dynamical Casimir effect. We speculate that a 
separated parametric oscillator for generating coherent GR microwaves could also be built. 

1 Introduction 

In this paper we discuss three apparently distinct phenomena: The gravitational Aharonov-Bohm ef- 
fect, the dynamical Casimir effect arising from parametric oscillations, and gravitational waves. The 
first of these phenomena is simply the gravitational version of the electromagnetic Aharonov-Bohm ef- 
fect. There has been recent interest in the gravitational version of the scalar Aharonov-Bohm effect 
IH. Here we will discuss the gravitational version of the vector Aharonov-Bohm effect. In the second 
phenomenon, i.e., the dynamical Casimir effect, we propose a possible experiment in which photons 
could be "pumped out of the vacuum" via a vibrating superconducting (SC) "membrane" considered as 
a parametric oscillator. Finally, we discuss gravitational waves. We speculate that a gravitational version 
of "pumping gravitons out of the vacuum" via parametric amplification, and above threshold, parametric 
oscillation, might be possible. The analog of a laser for gravitational waves could thus be constructed. 
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Figure 1: (Color online) Sketch of a gravitational Aharonov-Bohm effect. A "solenoid" with circulat- 
ing mass currents (in blue), produces "flux" (black dots) of a certain "gravito-magnetic" field (the Lense- 
Thirring field) in its interior (in yellow). In its exterior (in white), this field is zero. Nevertheless, an electron 
wave packet (in red), which is split at point A to go around the "solenoid" via paths 1 and 2, and then 
recombined at point B, will exhibit an Aharonov-Bohm fringe shift. 
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The concept that links these three phenomena together is the use of the DeWitt minimal coupling rule, 
whereby particles are coupled to both the electromagnetic and the gravitational vector potential. 



2 Gravitational Aharonov-Bohm effect 

A gravitational analog of the vector Aharonov-Bohm effect is depicted in Figure 1 ||T)||2|[3]. Aharonov- 
Bohm (AB) interference can occur when an incoming single-electron wavepacket is split at point A by 
means of a beam splitter into two partial waves traveling along paths 1 and 2, respectively, that go around 
the outside of a "solenoid" which contains circulating mass currents (indicated by the blue arrows). 

For instance, such mass currents in a cylindrical superconducting (SC) mass shell (indicated in black 
in Figure 1), could be produced by rotating the shell at a constant angular frequency around its cylindrical 
axis. The two partial waves could then be recombined into an outgoing single-electron wavepacket 
at point B by means of another beam splitter. Just like in the purely electromagnetic AB effect, the 
"flux" of a certain gravito-magnetic field (i.e., the Lense-Thirring field) would be confined entirely to the 
interior region (indicated in yellow in Figure 1) of the "solenoid," and would vanish at all points in its 
exterior (indicated in white). There results a quantum mechanical AB fringe shift due to the "flux" that 
is observable at point B, which cannot be explained classically. 

To understand the thought experiment pictured in Figure 1, we begin from Einstein's field equa- 
tions, which, in the limit of weak gravitational fields near slowly moving matter (i.e., in the vicinity of 
nonrelativistic masses), become the following set of four Maxwell-like equations |4|: 

V-E g = -^ (1) 
<9B„ 

V-B g =0 (3) 
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dE 

where the gravitational analog e g of the electric permittivity £o of free space is [5] 

1 

and where the gravitational analog /z„ of the magnetic permeability of free space is |6| 



VxB g = 4Hg[ -j g + E g ^r ) (4) 



e g = = 1.19 x 10 y SI units (5) 



AnG oi 
M„ = — r- = 9.31 x 10~ 27 SI units (6) 
c z 

Here G = 6.67 x 10~ n SI units is Newton's constant, and c = 3.00 x 10 8 m s in SI units is the vacuum 
speed of light. 

In the four Maxwell-like equations, ([TJl — (|4]), the electric-like field E g is the gravito-electric field, 
(i.e., the local acceleration g of a freely falling test particle), which could be produced by the mass density 
p g of nearby matter, via ([TJl. Likewise, the magnetic-like field B g is the gravito-magnetic field, which 
could be produced by the mass current density j g of nearby nonrelativistically moving matter, and also 
by the gravitational analog of the Maxwell displacement current density s g dE g /dt, via In the case 
of nearby stationary nonrelativistic mass currents, B s , can be identified with the Lense-Thirring field Q 
that is generated by these currents. 

A nonrelativistic test particle with a mass m moves in the presence of the weak fields E g and B ? in 
accordance with the Lorentz-like force law H 



-m—=m(E g +vxB g ) (7) 



where m is the mass of the test particle and v is its velocity (with v << c). 

To understand the experiment pictured in Figure 1, we only need the stationary version of Q, i.e., 
the gravitational analog of Ampere's law 

V x B, = -4ii g j g (8) 

where j g and B g do not depend on time. This implies the following gravitational analog of Ampere's 
circuital law: 

V^=-4M4)enc (9) 

where dl is a line element of an arbitrary closed curve C, and (I g ) enc is the mass current which is enclosed 
by C. Applying this Ampere's circuital law to the "solenoid" of Figure 1, which could be a uniformly 
rotating SC cylindrical mass shell, and using an appropriately chosen closed curve C, one concludes that 
the B g field in the interior to the "solenoid" is a uniform field pointing along the cylindrical axis, and that 
it has a constant magnitude 

B g = 4n g l' g (10) 

where I' g is the mass current per unit length of the "solenoid" flowing around the circumference of the 
rotating cylindrical mass shell. Furthermore, by another appropriate choice of the closed curve C, one 
concludes that everywhere outside the "solenoid," it is the case that 

B g =0 (11) 

i.e., that the Lense-Thirring field vanishes everywhere exterior to the "solenoid." This is analogous to the 
fact that the magnetic field vanishes at all points outside of an electromagnetic solenoid. Hence it follows 
from the Lorentz-like force law |7]i that although the electron in Figure 1 experiences a radial classical 
gravitational force due to the mass of the "solenoid," it could never have experienced any azimuthal 
classical gravitational force on its way from point A to point B via either path 1 or path 2 that could have 
caused the AB phase shift. 
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Put differently, if one thinks of the "solenoid" as a rotating SC cylindrical mass shell, the experiment 
has two independent parameters, namely, the linear mass density, and the angular velocity of the shell. 
That means one could shift the interference fringes by changing the angular velocity. Since the gravito- 
electric field from the mass of the shell does not depend on the angular velocity, a fringe shift will happen 
despite the fact that the classical force has not changed. Hence the AB fringe shift in the gravitational 
case could not have had a classical origin. 

Now from the Maxwell-like equation and from the vector identity 

V-(Vxh)=0 (12) 
it follows that it is always possible to express the magnetic-like field B ? as 

B s = Vxh (13) 



for some vector field h. The relationship ( 13 1 is formally identical to the relationship in electromagnetism 



between the magnetic field B and the electromagnetic vector potential A 

B = VxA (14) 

which follows from the Maxwell equation V • B = 0. Therefore we shall call h the "gravitational vector 
potential." 

In the gravitational case, just as in the electromagnetic case, the gravitational vector potential h 
possesses the gauge freedom 

h^h + V^ (15) 

where ji can be any arbitrary scalar function of position. This follows from the vector identity V x V/z = 
0, and is formally identical to the case of electromagnetism, in which the vector potential A possesses 
the gauge freedom 

A^A + VA (16) 

where A can be any arbitrary scalar function of position. Again, this follows from the vector identity 
V x VA = 0. 

Now the principle of local gauge invariance in nonrelativistic quantum mechanics states that the 
phase of the time-independent wavefunction *P (r) of any quantum system must always be able to be 
locally transformed without affecting the physics of the system. In other words, the transformation 1 8 1 

*P(r) ->¥(r)exp(z0(r)) (17) 

where the phase <j> (r) can be any arbitrary real scalar function of position r, can neither change the 
properties of the quantum system, nor the physical laws governing the system and its interactions with 
its environment. In particular, this local transformation of the phase of the wavefunction cannot change 
the probability distribution of the system, since 

\V (r) | 2 -> |¥ (r) exp (ty (r)) | 2 = (r) | 2 (18) 

and therefore the Born probability interpretation of the wavefunction is unaffected by this transformation. 

However, gradients of SP (r) will be changed by the introduction of an arbitrary scalar function 
(r), and therefore will alter the momentum of the system. If so, one could arbitrarily alter the physical 
laws governing the system, including altering the conservation of momentum of a particle in the usual 
exp (/p • r/h) plane wave state of an electron within a force-free region of space, where one knows that p 
must be a constant. This obviously cannot be the case. Therefore the principle of local gauge invariance 
necessitates the existence of some compensating vector field (or fields), such as the A and h fields in the 
DeWitt minimal coupling rule [9 | 

h h 

-V -> -V-qA-mh (19) 
i i 

P -> p-qA-mh (20) 
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where pop = \ V is the momentum operator, q is the charge, and m is the mass of the nonrelativistic quan- 
tum system under consideration. Here A and h are, respectively, the vector potentials for electromag- 
netism and for weak gravitation, which are being viewed here as being the requisite "compensating vector 
fields," whose existence is necessitated by the principle of local gauge invariance. Since the vector fields 
A and h have the gauge freedoms A (r) — >• A (r) + VA (r) and h (r) — > h (r) + V^i (r), these freedoms 
can then be used to compensate for the gauge freedom in the transformation *P (r) — > SP (r) exp (i<j> (r)), 
in just such a way that the quantum system can once again satisfy the principle of local gauge invariance. 

Thus invoking the DeWitt minimal coupling rule ( |20| , we demand that the nonrelativistic Hamil- 
tonian of any quantum system in the presence of A and h fields must always have the following form 
& 

1 ? 
H = — (p-qA-mhY+V (21) 
2m 

where V is the potential energy of the system. Here, in the present context of the SC quantum systems 
that we are interested in, such as that of the SC "solenoid" pictured in Figure 1, q = 2e is the charge of a 
Cooper pair, and m = 2m e is its mass. 

Although one can always arbitrarily choose a gauge locally so that both vector fields A and h are 
set identically equal to zero at each point exterior to the solenoid, nevertheless the fluxes interior to the 
solenoid 

4> = <j)A-dl (22) 

c 

4> g = jh-dl (23) 

where C is a closed curve enclosing the solenoid, cannot be arbitrarily set equal to zero, but must instead 
be gauge-invariant, nonzero, globally measurable quantities. Hence the fluxes <f> and 4> A , must be physical 
quantities. 

The gravitational Aharonov-Bohm effect depicted in Figure 1 is closely related to the time holonomy 
which arises from the off-diagonal time-space components of the metric tensor goi If X Oil . It can be shown 
ifm that this time holonomy At can be expressed as follows: 

At = —I—dx i (24) 
c J goo 
c 

where C is an arbitrary closed curve in space (such as the one enclosing the "solenoid" in Figure 1), 



and dx' is a spatial line element of this closed curve. In light of the time holonomy given by ( 24 1, it is 
impossible in general relativity to define a global time coordinate for an entire physical system, such as 
the topologically nontrivial superconductor in Figure 1 . 
In the weak field, slow matter approximation, in which 

got ~ hot (25) 

where h^ are the time-space components of the small-deviation metric h^ v from the Minkowski metric 
r\nv, and in which the time-time component of the metric can be approximated by 

Soo«-l (26) 



it follows that the time holonomy ( 24 1 becomes approximately 



At « - j hoidx' (27) 
c 



For electron waves traveling around a closed curve C enclosing a "solenoid" such as that in Figure 1, the 



time holonomy ( 27 i becomes the phase holonomy 



m c \ f 

A $ = ®Compton Af ~^ h --f hoidx ' ? (28) 

c 
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where G>c om pt n = rn e c 2 /h is the Compton frequency of the electron (T). The phase shift (28 1 , which is 
nonvanishing for the "solenoid" configuration of Figure 1 , is the gravitational AB phase shirt. It is closely 
related to Berry's phase [12], since both phases have a common origin in non-Euclidean geometry. 

Since physically counting the number of fringes in a shift of an Aharonov-Bohm interference pattern, 
such that in Figure 1, must yield the same result for all observers, independent of their state of motion 
under a restricted set of (Galilean) coordinate transformations, it is sufficient for the purposes of this 
paper to say that the flux is a Galilean invariant, and therefore a measurable, physical quantity. 



Moreover, the closed-path integral of h^dx' in ( 27 1 is an intrinsic time holonomy which cannot vanish 
due to an arbitrary gauge transformation. 

However, if one were to arbitrarily make the global gauge choice 

liQi = everywhere (29) 

as is done, for example, in the transverse-traceless gauge, then it follows that the phase holonomy must 
vanish identically, i.e., 

mc 2 1 / 

A(h ss <p liQidx' — by setting /to; = everywhere (30) 

n c J 

c 

and the gravitational Aharonov-Bohm phase shift predicted for the electron interference pattern in Figure 
1 would disappear. 

However, just as in the case of the electromagnetic Aharonov-Bohm effect where 

= H j> Ajdx' — by setting A; = everywhere (31) 



the results ( 30 1 and ( 3 1 1 are both unphysical whenever the closed curve C encloses either a solenoid with 
a nonvanishing electromagnetic flux 4> ^ 0, or a "solenoid" with a nonvanishing Lense-Thirring flux 
<t> g 7^ 0, since both of these fluxes are gauge-invariant, measurable, physical quantities that cannot be 



arbitrarily set equal to zero. Hence the transverse-traceless gauge choice ( 29 1 is unphysical in situations 
that involve the gravitational Aharonov-Bohm effect depicted in Figure 1 where <P g ^ 0, and also, by ex- 
tension, in time-varying situations that involve <P g (t) ^ 0, for example, in situations where gravitational 
radiation is interacting with superconducting systems. 



The usual Aharonov-Bohm phase shift follows from DeWitt's minimal coupling rule ( 20 1 when one 
sets <J> ? = 0, for then the phase shift arising from (20 1 in the configuration pictured in Figure 1, reduces 
down to the usual expression 

A0 = | jx-dl ^0 (32) 
c 

and we recover the standard form for the AB phase shift. However, if <P g ^ 0, then the total AB phase, 
upon integration over any closed curve C enclosing the "solenoid," becomes 

q f m f q<$> m4>„ 

A ^H? A - dl+ hV- dl =-h + ^ L (33) 

c 



The vector potential h in ( (33) can arise either from a Lense-Thirring field, or from rotations of a quantum 
system, such as from a rotating SC ring. Since experiments with rotating SC systems are much easier 
to perform than experiments involving Lense-Thirring fields, we shall confine our attention for now to 
these much easier experiments. 



One can check experimentally the expression given by (33 i for a rotating SC ring, in which case 
the single-valuedness of the macroscopic wavefunction of the Cooper pairs demands that A^oj = 2nn, 
where n is an integer. It follows from this that a magnetic field must be generated by rotating a super- 
conducting ring, i.e., that a "London moment" must accompany this rotation. Precision measurements 
of the London moment of a rotating SC ring thus provide a test for the correctness of the expression for 
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the total AB phase in ( |33j ). Cabrera and co-workers 1131 performed these measurements to 100 parts 
per million. Thus the formula in (|33]l for the total AB phase has been experimentally verified. In this 



way, the expression in (20 1 for the DeWitt minimal coupling rule has been experimentally tested to high 
precision. 

So far we have been considering only the case of stationary, time-independent, charge and mass 
currents, such as those in a SC magnet, or in a steadily rotating SC ring. These quantum currents can be 
the quantum mechanical sources of time-independent A and h fields that give rise to the AB effect. 



3 The dynamical Casimir effect via parametric oscillations 

In this section, we begin by discussing the case when the vector potentials A and h are time dependent, 
i.e. A(r,f) and h(r,f). From this we will propose a version of the dynamical Casimir effect for the 
electromagnetic vector potential A, in which photons are "pumped" out of the vacuum via parametric 
oscillations of a SC membrane. 

It is natural to extend the time-independent Hamiltonian pT) to the following time-dependent one 
El: 

H=^-(p-qA(r,t)-mh(r,t)) 2 +V (34) 

in which the fields A(r,f) and h(r,f) are to be first treated as classical fields, but the matter (e.g., the 
vibrating SC wire in Figure 2) is to be treated quantum mechanically, in the so-called "semi-classical 
approximation." 



Expanding the square in (34i, one obtains the following interaction Hamiltonian terms: 

#P-A = --P-A(r,f) (35) 
m 

which leads to the interaction of the quantum system with electromagnetic (EM) radiation, such as in the 
stimulated emission and absorption of EM waves by the quantum matter, and 

# p . h = -p-h(iy) (36) 

which leads to the interaction of the quantum system with gravitational (GR) radiation, such as in the 
stimulated emission and absorption of GR waves by the quantum matter, and 

H A . b = +qA(r,t)-h(r,t) (37) 

which leads to the interaction between EM and GR radiation fields mediated by the quantum system, 
such as in the transduction of EM waves into GR waves mediated by the quantum matter [ 15], and 

a 2 

ff A .A = +!^A(r,t)-A(r,f) (38) 

which leads to Landau-diamagnetism type of interactions of the quantum system with EM radiation, such 
as in the parametric amplification of EM waves by a strongly driven, i.e., "pumped," quantum system 
(see, for example, Figures 2 and 4), and 

m 

H h . h = +-h(r,t)-h(r,t) (39) 

which leads to gravitational Landau-diamagnetism type of interactions of the quantum system with GR 
radiation, such as in the parametric amplification of GR waves, again by a strongly driven, i.e., "pumped," 
quantum system (again, see, for example, Figures 2 and 4). 

All of the above interaction terms will be treated as small perturbations of the unperturbed Hamilto- 
nian 

P 2 

Ho = ^-+V (40) 
2m 
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and can thus be treated using standard perturbation theory. 

At a fully quantum mechanical level of description, both the matter and the radiation fields A and 
h would have to be quantized. The radiation fields could be quantized by invoking the commutation 
relations 

= 1 (41) 
M + ] = 1 (42) 

where a and a ' are, respectively, the annihilation and creation operators for a quantum of a single mode 
of the EM radiation field, and where b and are, respectively, the annihilation and creation operators 
for a quantum of a single mode of the GR radiation field. 

For now, let us focus solely on the interactions of quantized EM radiation with matter. The second 
quantized form for the EM vector potential operator Aop (r), when summed over all the modes of a 
cavity enumerated by the index K, is lfl6ll 




A op(r)=IW^—(fl, + 4)^(r) (43) 



where co K is the frequency of mode K, and S K (r) is one Cartesian component of the classical electric 
field distribution associated with this mode. Therefore, when there is a single dominant mode in the 
problem, the vector potential operator simplifies to the expression 

A op (r)c< (fl + flt) (44) 

where the mode index K and the proportionality constant have been suppressed 

The expansion ol 
contain the term ifTTll 

Kop oc a* a* + hermitian adjoint (45) 

The a* a* term of the operator Aop corresponds to the process of photon pair creation in the parametric 
amplification arising from the pumping action of some strong "pump" wave upon a quantum system. It 



The expansion of the square in the interaction Hamiltonian fl^ A x (a + a^) 2 in ( 38 1 will therefore 



can be shown that ( 45 1 has the form of an infinitesimal generator of a squeezed state of light [17]. 



Instead of enumerating all the possible resulting second quantized forms of the above interaction 



Hamiltonian terms, let us just focus on one such term, namely, //a a m ( 38 I, which is associated with 
parametric amplification, such as that in the setup depicted in Figure 2. This Figure represents an "opto- 
mechanical" parametric amplifier, which becomes, above threshold, a parametric oscillator, whose active 
element is the central vibrating SC wire (indicated in red), placed across the middle of an extremely high- 
Q SC microwave cavity. Here, instead of using optical cavities, as is usual in ongoing opto-mechanical 
experiments [ 18 1, we shall be using SC microwave cavities. The reason for this is that the quality factor 
for SC microwave cavities has already been demonstrated by Haroche and co-workers 1 19 1 to be on the 
order of 

Q ~ 10 10 (46) 

which can be much higher than that of typical optical cavities. 

The motion of the SC wire in the middle of the microwave cavity will modulate the "pump" mi- 
crowaves coming through the "IN" port so as to produce radiation at new sideband frequencies via 
the Doppler effect. The "seed" radiation initially in one of these sidebands, namely, the first "Stokes" 
sideband, can then become the exponentially amplified. Macroscopic, easily detectable radiation in the 
form of a strong Stokes wave emitted by the parametric amplifier can then leave the cavity through the 
"OUT" port. This kind of strong Stokes emission would be similar to the Stokes emission observed in 
the stimulated Raman effect in nonlinear optics 1201 . 

However, here we shall first analyze classically the parametric amplification process in Figure 2, in 
order to answer the following questions: What is the threshold for parametric oscillation in Figure 2? Is 
this experiment feasible to perform? The key concept that we shall use in this classical analysis is that of 
the work done by the moving wire, viewed as if the wire were a moving "piston" acting on some "seed" 
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Figure 2: (Color online) A parametric amplifier or oscillator, whose active element is the vibrating SC wire 
(in red) placed in the middle of a microwave SC cavity (in grey). The moving wire can be viewed as if it were 
an oscillating "semipermeable membrane," which does work upon some "seed" radiation initially present in 
the cavity, thus amplifying this radiation in a reciprocating, piston-like action. Photons incident upon this 
moving "membrane" experience a Doppler shift that changes their energy. Thus when pump microwaves 
enter through the left hole, an amplified signal (Stokes) wave will exit through the right hole. 



radiation initially present in the cavity. This "piston-like" action of the moving wire can also be viewed 
as if the wire were a partially reflecting "moving mirror," and will lead to the exponential amplification 
of the "seed" radiation at the Stokes frequency above the threshold of parametric oscillation. 

The reason for using a vibrating wire instead of a vibrating membrane is that the wire is one dimen- 
sional, whereas the membrane is two dimensional. The mass of a thin wire can be made much smaller 
than the mass of a thin membrane, and therefore a wire can be driven more easily into motion. 

However, here we shall model the vibrating wire in Figure 2 as if the wire were a vibrating "semi- 
permeable membrane," for which the pressure acting on the membrane can be converted into an easily 
calculable force. The justification for this "membrane" model is that the scattering cross section of a thin 
conducting wire, when placed symmetrically across the mouth of a waveguide, can be comparable in 
size to the cross-sectional area of the waveguide, because the wire tends to "short out" the electric field 
of the TE mode of a waveguide. Thus the reflection coefficient of the wire placed across the middle of a 
cavity (as in Figure 2), can be made quite high. There results a splitting of a microwave cavity mode into 
a spectral doublet as illustrated in Figure 3 |F2Tlll22l . due to the presence of the central wire. The splitting 
frequency of the doublet can typically be on the order of 1 GHz for a microwave mode frequency of 
around 10 GHz. 

We have observed the splitting of a microwave cavity resonance into a spectral doublet. An RF 
cylindrical copper cavity of length L = 1 .284" and diameter D = 1 .02" with two parallel conducting end 
plates was constructed to support a TE1 12 mode at 11.42 GHz 11231 . (The mode indices l,m,n are chosen 
to correspond to the number of half wavelengths along their respective axes; angular, radial, and axial 
respectively.) The input coupler, a short straight wire placed perpendicular to the axial direction, was 
placed on the cylinder at approximately a quarter wavelength from one end plate. The output coupler is 
a small loop placed at the other end plate of the cavity. 

An off-diagonal scattering matrix element S21 transmission measurement, preformed with a HP 
8720C network analyzer, shows the resonance of the TE112 mode at approximately 11.50 GHz; see 
figure |4|a). A splitting is observed by placing a 22AWG copper wire at the midpoint of the cavity, 
perpendicular to the axial direction and parallel to the input coupler. The splitting is approximately 400 
MHz; see figuregb). 
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Figure 3: The excitation spectrum of a microwave cavity with a wire placed at its center (see Figure 2). The 
spectrum consists of a doublet of resonances at the "signal" mode at the "Stokes" frequency (Og and the 
"pump" mode at the "pump" frequency (O p . The difference between "pump" and "Stokes" frequencies is 
resonant with the frequency Q. of the vibrating wire (i.e., co p — (Og =£1). The "anti-Stokes" frequency (Oj^ 
is off resonance with respect to the doublet, and hence is suppressed. (Cf. the stimulated Raman effect in 
EQl). 



In Figure 3, the pump frequency 03 p of the parametric amplifier is assumed to be tuned to coincide 
with the upper member of the spectral doublet, and the signal frequency 0)g is assumed to be tuned to 
coincide with the lower member of this doublet, which we shall call the "Stokes frequency," in analogy 
with the stimulated Raman effect [20|. The idler frequency ft),, i.e., the frequency of the mechanical 
motion of the central wire in Figure 2, is the beat frequency Q. = co p — (0§ between the pump and signal 
frequencies. Note that the parasitic, Doppler upshifted "anti-Stokes frequency" C0j^ is automatically 
suppressed by this spectral doublet. 

To calculate the force acting on the centtal membrane (as a model of the force acting on the centtal 
wire) in the middle of the microwave cavity of Figure 2, we begin from the Maxwell stress tensor [24 1 1 25 1 

1 „ 9 \ If 1, , 



Tij = e \EtEj - -8 i} E 2 \ + - [BiBj - -8^) (47) 

Now starting from the electromagnetic force exerted on charges and currents given by the Lorentz 
force law 

F = c/(E + vxB) (48) 

it can be shown that there results the following relationship between the force F and the total Maxwell 
stress tensor 7y and the Poynting vector S [24||25 1: 

r J r 

(49) 



(F), = j ■ (da) J - eoMo^ J (S),</V (where i = 1,2,3) 



S(V) v 

where 7},- is the stress tensor evaluated at da, an infinitesimal area element of an arbitrary surface S (V) 
which encloses the volume V, dV is an infinitesimal volume element of the matter inside the volume V 
enclosed by the surface S (V), and S = E x H is the Poynting vector evaluated at dV inside V. 

Now the tangential electric field must vanish at the boundary of any conductor. Hence, for all trans- 
verse electric modes of the microwave cavity pictured in Figure 2, if one chooses the surface S (V) to be 
that of a small pillbox straddling a patch of the surface of an equivalent SC membrane, the contribution 



to the force ( 49 1 from the Poynting vector term evaluated at the pillbox enclosing the patch of the surface, 



must vanish. In the case of transverse magnetic modes of the cavity, the Poynting vector S = E x H does 
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Figure 4: Splitting of the TE112 mode in cavity with a bisecting copper wire at its midpoint perpendicular 
to the axial direction and parallel to the input coupler. The splitting is on the order of 400 MHz. S21 
transmission measurements (a) of a copper cavity with TE112 resonant frequency at 11.5 GHz and (b) 
splitting due to the placement of copper wire placed at the center. The vertical axes use the same arbitrary 
power reference in the conversion from logarithmic to linear scale. 



not vanish at the surface, since there will be a longitudinal component of the electric field at the surface 
along with a tangential component of the magnetic field. Hence there will arise a tangential component 
of S at the surface of the SC membrane, but this S cannot contribute to any normal force acting on the 
conducting surface. 

Therefore the only contribution to the normal force acting on the equivalent SC membrane arises 




solely from the Maxwell stress tensor term of (49 1, which, for the case of transverse electric modes 
evaluated at the membrane, reduces down to 

(Tij) - — ( +B 2 y | (50) 

because the magnetic field of the transverse electric mode (e.g., the TE112 mode), whose electric field 
is pointing in the x direction (in the Cartesian coordinate system of Figure 2), will be pointing in the y 
direction at the surface of the equivalent membrane l26ll . 

If one therefore replaces the SC wire by an equivalent SC membrane, then the diagonal terms of ( |50| 
can be interpreted as a "field pressure" acting on the membrane with a maximum amplitude of 

^ax-^(B 2 ) max = M max (51) 

where (B 2 ) may \s the maximum of the square of the magnetic field, and («B) max is the maximum magnetic 
energy density, evaluated at the surface of the membrane. 

Because the stress tensor depends quadratically on the field, there results a pressure being exerted 
on the membrane at a beat note frequency due to the beating between the fields at pump frequency co p 
and the fields at the Stokes frequency a>§ in the spectral doublet of Figure 4. This beat note can drive 
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the membrane (or the wire) at the beat frequency £2 — co p — 0)§, i.e., at the splitting of the upper and 
lower members of the doublet. The force in the z direction acting on the membrane (or wire) at the beat 
frequency £2 will therefore have the form 

F a = — 3§ p 3§~ exp (-i ((Op- (O s )t)-£/ e te + c.c. 
Mo 

oc exp(-/£2f) +c.c. (52) 

where £^ e ff is an effective area of the membrane (or, equivalently, the effective scattering cross-section 
of the wire), and where 

Bp — ^ p exp(-/oy) +c.c. (53) 

is the pump waveform, with SS p being the complex amplitude for the pump magnetic field waveform, 
and where 

B s = ^ s exp (-j'OHj* ) + c.c. (54) 

is the Stokes waveform of some small amount of "seed" radiation already present inside the cavity, 
with 8§v being the complex amplitude for the Stokes magnetic field waveform. (Note that such "seed" 
radiation could in principle be vacuum fluctuations of the EM field inside the cavity.) In the expression 
( 52 1 for the force Fa at the beat frequency £2, we have assumed that the pump wave ( 53 1 is always much 
stronger that the "seed" Stokes wave (p4|i, i.e., \3§p\ » l^c . 

Since the driving force Fa at the beat frequency £2 can be made resonant with the acoustical resonance 
frequency of the membrane, we shall model the resulting motion of the membrane as that of a simple 
harmonic oscillator with a resonance frequency of £2. Using Newton's equation of motion for a damped 
simple harmonic oscillator moving in the z direction, viz., 

m(z + yz + £l 2 z) =F n (55) 

where m is the mass of the membrane, and y is its damping coefficient, and using an Ansatz of the form 

z = Zmax exp (-**£&)+ c.c. (56) 

for the displacement of the membrane in the z direction, one finds that its maximum, on-resonancce 
complex displacement amplitude is 

Zmax = i J^eff ( 57 ) 

The velocity of the membrane in the z direction will then have the form 

v = v max exp (— iQl) + c.c. (58) 
where the complex velocity amplitude v ma x = — i&Zmax becomes, on resonance, 

Vmax = -^ e ff (59) 

Homy cn 

There results a Doppler effect arising from the velocity of the membrane moving in the z direction in 
Figure 2, giving rise to upper and lower Doppler sidebands. However, only the lower Doppler sideband 
will be excited, since only the lower sideband will be resonant with lower member of the doublet in 
Figure 4. 

The maximum, on-resonance time-averaged power (P) being delivered from the pump wave into 
the simple-harmonic membrane motion at the beat frequency £2, and therefore into the lower Dopper 
sideband, i.e., into the Stokes wave ( f54] i, is 

1 



= (Fa ■ v) raax - 2 \^ p \ z |% \ 2 -2— < ff (60) 
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Now invoking the conservation of energy, we find that, if we for the moment neglect all losses, the 
Stokes wave will be amplified by this power transfer, such that the power (P m ax) being tranferred into 
the Stokes wave must equal the rate of growth of the time-averaged Stokes energy inside the cavity 

V^s) = ^ l%| 2 ^ff = ^ l%| 2 ^eff^eff (61) 

where -S? e ff is an effective length of the cavity (i.e., V e ff = ^ffJz^eff is an effective volume of the cavity). 
In other words, from energy conservation it follows that 

(Pmax)-^(f/ S ) ( 62 ) 

Substituting in from ( |60) and (61) , one infers that 

2 K l*S I' ^T/eff = ! 1% I' ^eff) - ■* 1% I 2 ^eff) («) 

where Kg is the exponential gain coefficient for parametric amplification of the Stokes wave. Thus we 
arrive at an exponential-growth ODE for the energy stored U§ in the cavity at the Stokes frequency 

J t (Us) = Ks(Us) (64) 



Solving for the gain coefficient K§ from ( 63 1, we conclude that 

_ 2\,% p \ 2 £/ eS 

Thus the gain of the Stokes wave is directly proportional to the pump power stored in the cavity, just like 
in the stimulated Raman effect [20|. 

Like in a laser, the threshold for oscillation [27 1 1 28 1 occurs when 

Gain = Loss (66) 

Above threshold, i.e., when the gain exceeds the loss, macroscopic amounts of coherent radiation can 
be produced inside the cavity by the exponential amplification of the "seed" radiation, i.e., of vacuum 
fluctuations. Thus even if one were to start off only with vacuum fluctuations as the "seed," one can 
produce macroscopic amounts of coherent radiation by the stimulated emission of radiation, just like in 
a laser. In principle, this should also apply to GR radiation, as well as to EM radiation. Thus generators 
of microwave frequency GR radiation should in principle be possible to construct, as well as amplifiers 
and detectors for this kind of radiation. The only remaining question is whether such devices are feasible 
in practice. 

The loss of the Stokes wave (i.e., the "signal") from the cavity depicted in Figure 2 can result from 
emission of radiation through an outcoupling hole into the environment, or from remnant ohmic losses 
in the components of the cavity. We shall call the resulting quality factor of the cavity at the Stokes 
frequency under working conditions (i.e., including internal losses and the outcoupling into the environ- 
ment) the "loaded Q", which is defined as follows: 

2s =<0 S T S (6 7 ) 

where "S" stands for "Stokes wave" whose frequency is (Q§, and whose stored energy inside the cavity 
decays away with a time scale Tqj (the so-called "cavity ring-down time") after the pump has been turned 
off. 

Moreover, the loss coefficient y in the motion of the simple harmonic oscillator leads a decay time 
T n = 1 1 J f° r the energy stored in the oscillator. This leads to a mechanical oscillator quality factor Qq, 
which is defined as follows: 

Qa = £}r n (68) 
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It follows from the gain-equals-loss condition d66b that at threshold 



t K \ = 2 I I threshold ^eff = 2 1 I threshold Qn-^eff = 1 = ®S rfi ™ 

I S; threshold Homy^g Lk>m£l3? eS t s g s V ; 

Since the time-averaged stored energy in the pump wave inside the cavity depicted in Figure 2 is 

(^) = ^|^ P | 2 <ff^eff ( 7 0) 



we conclude from ( 69 1 that the threshold pump power needed for parametric oscillation is 



(^threshold = 2 Q s Q n (71) 

This result is to be compared with the threshold pump power needed for parametric oscillation for ex- 
citing an elastic mode of a mirror of a Fabry-Perot resonator obtained by Braginsky and co-workers 
ED 

in \ - 1 m(0 ^ Ll n2 ) 

\ u Pl Braginsky ~ 2 Q t Q s 

where m is the mass of the mirror, CO s is the frequency of the elastic mode, L is the length of the Fabry- 
Perot resonator, Qj is the quality factor of a down-shifted "idler" optical mode of the resonator, and Q s 
is the quality factor of the elastic mode. By inspection of ( f7T| > and (72 1, we see that these thresholds are 
quite similar. 

However, the electrodynamic Q factor of SC microwave cavities is typically on the order of 10 10 lfT9l . 
whereas the typical mechanical Q factor for the best opto-mechanical oscillators, which are composed 
of non-SC materials in the ongoing opto-mechanical experiments, is at most on the order of 10 5 |18|. 
Therefore the question naturally arises whether it is possible to replace these low-g, non-SC mechanical 
oscillators, with high-g SC mechanical oscillators, in which their mechanical Q can approach the typical 
electrodynamic Q ~ 10 10 of SC microwave cavities. 

One possible answer to this question is the "triple SC microwave Fabry-Perot resonator" shown in 
Figure 4 ET1 . in which a charged SC membrane is extremely tightly coupled via its electrostatic charge 
to the longitudinal electric fields of a transverse magnetic mode of a high Q SC microwave cavity. 
The idea here is that when the charge on the SC membrane is sufficiently large, then the mechanical 
dynamics of the membrane will be "slaved" to follow closely the electromagnetic dynamics of the high- 
Q microwave SC mode. Calculations ll29l show that one only needs a charge of pico-Coulombs for 
this to happen. Note that as a result of the "slaved" dynamics, the SC membrane will be moving at 
microwave, and not at acoustical, frequencies. This means that the motion of the membrane will be 
essentially that of a "free" mass, which is being driven solely by Maxwell's stress tensor. Therefore this 
microwave-frequency motion will be independent of the elastic and dissipative mechanical properties of 
the membrane. 

Another important feature of the configuration shown in Figure 4 is that the signal and the idler 
waves are spatially separated into two disjoint, high Q SC cavities, which are separated from each other 
by a common, vibrating SC membrane. It turns out that this leads to two separate Q factors in its 
denominator of the threshold for parametric oscillation, which arises due to this separation. We shall 
therefore call the parametric oscillator configuration of Figure 4 a "separated parametric oscillator," in 
contrast to that in Figure 2, which we shall call an "unseparated parametric oscillator." The threshold for 
the separated parametric oscillator of Figure 4 will turn out to be at least a factor of 10 5 lower than that 
of the unseparated parametric oscillator of Figure 2. 

We start the analysis of the separated parametric amplifier depicted in Figure 4 by examining the 
work done by the "pump" wave on the moving SC "membrane," when it produces a displacement by an 
amount Az of the membrane to the left along the axis of the "double" Fabry-Perot resonator on the right 
side of the membrane. The work done during this displacement is 

AW=( — B 2 ) • j/ eff Az (73) 
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Figure 5: (Color online) A "triple" SC microwave Fabry-Perot cavity [21] consists of a "single" SC cavity 
separated by a vibrating SC "membrane" (in red) from a "double" SC cavity with a fixed SC iris (in black) 
at its center. This separating membrane is totally impermeable to all microwaves. The membrane is electro- 
statically charged on it left surface with a charge +q, and the left SC mirror is charged on its right surface 
with a charge — q. "Pump" microwaves (in blue) enter into the system through the right hole, and "signal" 
waves leave the system through the left hole, but "idler" waves leave the system through the right hole. 



where is the effective hemiconfocal spot size at the membrane, Az is the displacement of this 
membrane, and 

ub = ^-B 2 (74) 
is the energy density of the magnetic field evaluated at the right surface of the "membrane," which is the 



pressure arising from the Maxwell stress tensor (50 1, i.e., a pressure being exerted upon the membrane 



that can cause a change of the volume inside the "double" Fabry-Perot resonator on the right side of the 
membrane 

AV = £/ eS Az (75) 

where is the effective area of the membrane, which is determined by the hemiconfocal spot size of 
the mode on the right side of the membrane. For simplicity, we shall assume here that the hemiconfocal 
spot size of the mode on the left side of the membrane also has the same ■ 

The instantaneous mechanical work AW in ( |73] l done by the "pump" upon the fields of the resonator 
can be rewritten in the form 

AW = PAV (76) 
where the instantaneous pressure P on the "membrane" is 

P=-^B 2 (77) 

which is equal to the instantaneous energy density ug in d74b evaluated at the right surface of the mem- 
brane. It is clear from the expression for the work in (|76| that AW can be interpreted as if it were the 
work being done by a moving piston acting on a thermodynamic system, here, the radiation fields inside 
a cavity. 

Now we shall presently see that if energy were to be continually supplied to the "double Fabry- 
Perot" resonator on the right side of the membrane by some continuous-wave, external microwave pump 
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waveform oscillating at a frequency co p (i.e., the "pump" frequency) entering through the right hole 
of Figure 4, then the exponential amplification of some seed "signal" waveform at a frequency of co s 
within the "single" Fabry-Perot on the left side of the membrane, simultaneously with the exponential 
amplification of some seed "idler" waveform at a frequency of 0), within the "double" Fabry-Perot on the 
right side of the membrane, can occur. This amplification effect can arise from the mutual reinforcement 
of the signal and idler waves at the expense of the pump wave, in which the pump wave beats with 
the idler wave via the Maxwell stress tensor to produce more of the signal wave, and the signal wave 
modulates the pump wave via the Doppler effect to produce more of the idler wave, etc. The mutual 
reinforcement of the two "seed" waves will lead to an instability above a certain threshold, i.e., to the 
parametric oscillation of both the signal and idler waves that produces macroscopic amounts of both 
kinds of waves, which then leave the system in opposite directions via the left hole and the right hole of 
Figure 4, respectively, just like in a laser. 

For parametric amplification to occur, the frequency-matching condition 

(O p = co s + (Oj (78) 

must be satisfied. The meaning of the relationship can be most easily seen by multiplying it by the 
Planck's constant Ti so that one obtains the relationship 

h(o p = h(o s + hcoj (79) 

In other words, in the parametric amplification process, one signal photon hco s is simultaneously created 
along with one idler photon hco, at the expense of one pump photon hco p , which is annihilated during this 
"photon pair-creation process." In this process, an entangled pair of signal and idler photons, with the 
signal photon appearing on the left side, and the idler photon appearing on the right side of the membrane, 
respectively, will be produced in a correlated emission event inside the "triple" Fabry-Perot resonator 
depicted in Figure 4. This photon pair-creation process is described by the interaction Hamiltonian 

i/ mt oc a p a\a] + hermitian adjoint (80) 

which is a generator of a two-mode squeezed state 1 16 1. 

Let the microwave pump magnetic field just outside of the right surface of the membrane have the 
form 

Bp — & p exp(—i(Opt) +c.c. (81) 

where the pump magnetic field vector points transversely to the membrane immediately outside of its 
right surface, with SS p being the complex amplitude of the pump magnetic field. 

Similarly, let the "seed" idler magnetic field just outside of the right surface of the membrane have 
the form 

Bj = ^exp (-i(Ojt) + c.c. (82) 

which is a vector parallel to the magnetic field vector of the pump wave immediately outside of its right 
surface, with SSi being the complex amplitude of the idler magnetic field. We shall assume that the pump 
is tuned to be on resonance with the upper member of the spectral doublet of the "double" Fabry-Perot 
(see Figure 3), and that the idler is tuned to be on resonance with the lower member of this doublet. 
However, the off-resonance, parasitic "anti-Stokes" (i.e., the Dopper up-shifted) frequency component 
arising from the motion of the membrane will be suppressed, and hence neglected. 

To calculate the coefficient of parametric amplification, let us assume that the pump wave is much 
stronger than both a very weak "seed" idler wave and a very weak "seed" signal wave, so that LSt?J >> 



| S§i | and | SS p | > > \SS S \. It follows from ( 78 1 and ( |82| that the square of the total magnetic field evaluated 
at the right surface of the "membrane" will have the form 

B 2 = (Bp+Bjf = exp (-/ay) +^,exp(-/a),f) +c.c.) 2 (83) 
If we define the "beat frequency" as 

£2 = (Op - a>i (84) 
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then we see that there will arise cross terms in the square of the magnetic field ( |83] l which will contain 
terms that vary at the beat frequency £2, viz., 

(B 2 ) Q = & p exp(-ieo p t)&* exp(+ieOit)+c.c. 

= @p^texp(-iQt) + c.c. (85) 

Therefore there will exist a pressure being exerted on the membrane that varies at the beat frequency £2 
of the form 

{P)a = 2^0 = i (^^* ex P(- /i2f ) +cc.) (86) 

If we define the complex pressure amplitude &q as follows: 

= ^-Stp&i (87) 

then the pressure exerted on the membrane which varies at the beat frequency £2 will have the form 

(P) a = & n exp (-Kit) + c.c. (88) 

But the beat frequency £2 will be assumed to be tuned into resonance with the signal frequency, i.e., 

£2 = (Dp - to, = co s (89) 

so that the membrane can be driven at the resonance frequency co s of the "single" Fabry-Perot resonator 
to the left of the membrane. Thus power from the right side of the membrane can be fed resonantly by 
the motion of the membrane into the signal "seed" waveform on the left side. 

Due to the presence of the electrostatic charge +q on the left surface of the membrane, and the 
relationship 

F t (t) = qE z (t) (90) 

where F z (t) is the Coulomb force on the membrane exerted on the charge +q by the longitudinal electric 
field E z (t) of a transverse magnetic mode of the "single" Fabry-Perot resonator on the left side of the 
membrane [21 ], it follows that the motion of the membrane in the longitudinal z direction will be "slaved" 
through E z (?) to the dynamics of this transverse magnetic mode ll2~Tl Appendix B, where it was shown 
that one only needs q w 20 pC for the Coulomb force to dominate the dynamics of the membrane]. This 
is due to the tight coupling between Az(t) and E z (t) which arises from the electrostatic charge +q l29ll . 
Hence let us introduce an Ansatz that the displacement of the membrane has the form 

(Az(f)) n = £a(0 ex P(-' £i 0+ c - c - (91) 

where (t) is some slowly-varying complex displacement amplitude of the membrane, which is the 
slowly-varying envelope of the fast beat frequency phase factor exp (— iQl), 
Hence the velocity of the membrane will have the form 

(v)n = ^^ a «vn(0«ap(-in») + cc. (92) 
where the slowly-varying complex velocity amplitude of the moving membrane is 

v n (?) w -j£2e Q (r) (93) 



within the "slowly varying envelope approximation" [30], and the acceleration of the membrane will 
have the form 

^ n= (S) ~ a "( f ) ex P(" /£2f )+c-c (94) 
where the slowly-varying complex acceleration amplitude of the moving membrane is 

an(t)^-£l 2 e a (t) (95) 
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also within the slowly varying envelope approximation. 

The force due to the pressure ([86| being exerted on the membrane will have the form 



(F) a = & a (t) exp (-iQt) + c.c. 
where the slowly-varying complex force amplitude acting on the membrane is 



(96) 



(97) 



where is the effective area of the membrane, and where, in the "undepleted pump" approximation 
||20l , we have assumed that the pump amplitude 38 p is independent of time, but that the idler amplitude 
33* (?) can be a slowly varying function of time due to its amplification. 

Then the time-averaged mechanical power fed into the membrane's motion from the radiation, and 
hence into the signal wave of the "single" Fabry-Perot cavity of Figure 4, using ( 93 1 and ( 97 1, is 

(F-v) 



dW\ 

dt I signal 



c.c.) = (J? Q • i£lE a + c.c. 

1 



-2Im( i r Q Q£&) = -2Im 



Mo 



(98) 



Let the complex amplitudes of the pump, idler, and signal waveforms have the complex polar forms 



en 



|^/|exp(ifc) 
|e n |exp(/0. s ) 



(99) 
(100) 
(101) 



By inspection of ( 98 1, we see that the maximum power transfer from the radiation fields into the mem- 
brane's motion occurs when the phases of pump, signal, and idler waveforms (i.e., ( 99 1, ( 100 1, and ( 101 1) 
are adjusted so as to satisfy the condition 



<Pi ~ 0* = 



whereupon the maximum mechanical power fed into the membrane becomes 

\m p \\m(t)\\ea(t)\ 



/dW\ 
\~dF/ 



2^ eff n 



max, signal 



Mo 



(102) 



(103) 



Neglecting for the moment all dissipative losses, the kinetic energy of the membrane must grow due 
to this positive mechanical power being fed into it. Hence, invoking the principle of the conservation of 
energy, we get the equation 



/dW\ 

\ dt / max ^ s ig na i at 



(104) 



where we have used (93 1 and the fact that (v 2 ) = 2 1 Va| 2 - Therefore 
(m£l 2 \e n \ 2 )=2mtf\en ,dl£a 



d 
dt 



2^ eff n 



dt 



Mo 



(105) 



We thus arrive at an ODE for the rate of growth of the magnitude \e&\ of the displacement of the mem- 
brane 



d\eq\ 
dt 



■^eff 



(106) 
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Next, we shall obtain a similar ODE for the rate of growth of the magnitude |^, | of the idler wave. 
We start from the motional EMF created by the motion of the vibrating SC membrane, which leads to 
the generation of the motional electric field 

E = vxB (107) 

This relationship implies that the sinusoidal, back-and-forth motion of the mirror at a frequency Q. will 
be modulating the magnetic field oscillating at the pump frequency co p , such that an idler electric field at 
the surface of the mirror oscillating at the idler frequency 0), will be generated, i.e., 

(E) ( = (v) n x(B) p (108) 

This is a manifestation of the Doppler effect, in which the sinusoidal motion of the mirror will produce 
Stokes and anti-Stokes sidebands around the pump frequency co p . However, due to the doublet spectrum 
depicted in Figure 3 in the "double" Fabry-Perot resonator, only the down-shifted, first-order Stokes 
sideband will be resonant with the resonator. Therefore we shall neglect henceforth the anti-Stokes 
sideband, and all the other higher order Doppler sidebands. 

Note that when the membrane is moving towards the left in Figure 4, which is the direction in which 
the magnetic pressure due to the pump wave is pushing, this pressure will deliver power into the motion 
of the moving membrane, and simultaneously, will deliver power into the red-shifted, first-order Doppler 



sideband, i.e., the idler wave, via the relationship ( 108 I. This will lead to parametric amplification of the 
membrane's motion. 

Using the Cartesian coordinate system shown in Figure 4, we shall assume that the instantaneous 
velocity of the mirror is pointing in the — z direction, and that the instantaneous pump magnetic field 
vector Bp is pointing in +x direction, so that the instantaneous motional E field will be pointing in the 



+y direction. Thus, in terms of the complex amplitudes, ( 108 i reduces down to 



Si = v* a 3§ p (109) 

Now the time-averaged power delivered into the idler wave by the motional E field acting on the 
current density j induced by the idler wave at the right surface of the SC mirror, is 

dW\ f 

— ) = / <J • E) dV = 014) ^eff 5 + c - c - (HO) 
al I idler J 

where j* is the complex conjugate of the supercurrent density amplitude flowing on the surface of the 
mirror at the idler frequency (Oj, <£{ is the complex motional electric field amplitude at ©,, JZ^ff is the 
effective focal area of the membrane, and 8 is the London penetration depth, within which the supercur- 



rents j* will be flowing near the surface of the membrane. To calculate ft in ( HOi, we use Ampere's 



circuital law and a rectangular loop straddling the surface of the SC mirror, to get 

^*=fioj*8 (111) 



where 8 is London's penetration depth. Solving for j* from (JTTTJ, and substituting it into (110 1 using 
( |109| l, we get 



/^\ = ^(v^^;x ff +c.c. 

\ dt / idler Mo 



— ((«Qe£)^,^f)^ eff -|-c.c. (112) 



where in the last step we used ( 93 i for the complex conjugate of the complex velocity amplitude of 
the membrane. 



To maximize the power transferred to the idler, we again choose the phase condition ( 102 1 between 
the pump, idler and signal complex amplitudes, and find 

I dW\ 2 

( — ) =+-n|e Q | \m p l^-Kff (H3) 

\ at / maXi idler IM) 
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Assuming the absence of all dissipation, and invoking once again the principle of the conservation 
of energy, but this time for the idler wave, we obtain 



/dW\ 

\ at I max idler 



•^eff-^eff 



Mo 



-G|en| 



y ;|^eff 



(114) 



where Jz? e ff is the effective length of the "double" Fabry-Perot resonator. We thus arrive at an ODE for 
the rate of growth of the idler wave 



><' = 

where the constant of proportionality K% is 



eff 



\ea\\& P \ =K 2 \e Q \ 



& i 

K2 = -&-\9 p \ 

-^eff 

Note that this is of the same form as the ODE for the rate of growth of the signal wave 
earlier in ( |106| l, viz., 



where the constant of proportionality K\ is 



Ki 



eff 



(115) 

(116) 
obtained 

(117) 
(118) 



This implies that there exists a mutual enhancement of the signal and idler waves that leads to their 
exponential growth, i.e., to the parametric amplification of both waves. To see this, let us rewrite the two 



equations ( 1 15 1 and ( 1 17 1 in the following 2x2 matrix form: 



d 
dt 



Ki 

K 2 



= A 



(119) 



where A is the eigenvalue of the 2x2 matrix, viz.. 



*eff I 



\ jUo^eff 



The solution of ( |1 19| is 



leal 



leal 



exp (At) 



(120) 



(121) 



f=0 



The meaning of the positive root for A is that it represents the rate of exponential growth of the am- 
plitudes of the coupled signal and idler waves, when the phase condition ( 102 1, <p p — fa — fa = —Tt/2, 



for maximum power delivery to these coupled waves, is satisfied, whereas the negative root for A is 
that it represents the rate of exponential decay of the amplitudes of the coupled signal and idler waves, 
when the anti -phase condition, fa., —fa — <j) s = +7t/2, for maximum power extraction from these coupled 
waves, is satisfied. Whether one gets exponential growth or exponential decay of the waves thus depends 
on the choices of the initial phases of the pump, signal, and idler waves. This kind of phase-dependent 
amplification is the signature of the production of a squeezed state of the vacuum. 

Next, let us introduce a dissipative loss phenomenologically into the ODE for the idler \\\5\ as 
follows: 

^m\--\@i\=K2\en\ (122) 
dt Tj 
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where T; is the "cavity ring-down time" for the energy stored in the idler cavity mode on the right side 
of the membrane after the pump wave has been suddenly shut off, and, similary, into the ODE for the 
signal dTTTb as follows: 

-d<to|--|£h|=*l|#i| ( 123 > 
at T s 

where T s is the "cavity ring-down time" for the energy stored in the signal cavity mode on the left side 
of the membrane after the pump wave has been suddenly shut off. The relationships between the cavity 
ring-down times T, and T v of the two cavity modes, and their loaded quality factors Qi and Q s , are 

Qi = (Oi-d (124) 
Qs = 03sT s (125) 

At the threshold of parametric oscillation, there is a balance between gain and loss such that there 
arises a steady-state situation in which 

^1^1 = ^1^1=0 ( 126 > 
Therefore, at threshold, the two ODE's ( 122) and ( 123 1 reduce down to the two algebraic equations 

--\&t\ = K 2 \e a \ (127) 

-~M = Ki\m\ (128) 

*s 

Multiplying the left sides and the right sides of these two equations together, we get 

— =K 2 K, (129) 



By using the relationships ( |1 16| l, ((TTBj), ( |124| >, and ( |125| l, we get from ( | 129 1 the threshold condition 

(130) 



QiQs IM)m^f eff Hom^ {{ 
Since the time-averaged stored energy stored in the pump cavity mode is 

1 /,2W 1 ^ ,2, 



^ = 2^^ V ^ = ^ p \ Veff (131) 



we arrive from ( |130| ) at the conclusion that the threshold condition is 



\ U p /threshold = ( 1 32 ^ 



This is to be compared with Braginski's threshold condition (f72]> [27 1 



, x Braginski _ 1 mco^L 2 

\ u p /threshold ~ 2 Q t Q s K ' 

The above two expressions agree as to an order-of-magnitude estimate for the threshold of parametric 
oscillation for the "triple" Fabry-Perot resonator configuration of Figure 4. 

The required threshold input pump power (Pp)threshold ^ or P arameu "i c oscillation due to pump 
microwaves entering in through the right hole of the "triple" cavity configuration of Figure 4, can be 
found via the steady-state condition 

( P p)threshold = T~ ( U p) threshold < 134 ) 
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where the quality factor for the pump cavity mode Q p is related to the pump cavity ring-down time x p by 

Q P = co p t p (135) 



Finally, putting this together with ( 132 1, we conclude that for parametric oscillation to occur in the 
configuration of Figure 4, we need to inject a microwave pump power into the "triple" Fabry-Perot 
cavity the minimum amount of 



/threshold — 
Numerically, if we assume that ITJT 



(^threshold - q QiQs 6 (136) 



m — 2 mg (137) 

(O p =2%x 20 GHz (138) 

®i « (O s « 2k x 10 GHz (139) 

jSf eff sa Aj /2 « k s /2 « 3 cm (140) 

e P « a « a ~ io 10 (wd 

then we conclude that we would require a microwave pump power at a frequency of 20 GHz to be 
injected through the right hole of the "triple" Fabry-Perot cavity of at least 

< P p)threshold ~ 4 microwatts (142) 

Above this minimum power level, there would result a parametric oscillation effect in which a macro- 
scopic amount of signal and idler microwaves centered around 10 GHz, with powers on the order of 
microwatts (i.e., with powers comparable to the pump threshold power), would be emitted in opposite 
directions through the left and the right holes, respectively, of the "triple" Fabry-Perot cavity, like in a 
laser. If the Q is lowered by opening the outcoupling holes, the threshold will go up, but so will the 
output power of the parametric oscillator. For example, by lowering all the Q's to 10 9 instead of 10 10 , 
the threshold will be increased to 400 microwatts, but the output power of the dynamical Casimir effect 
will also increase by a few milliwatts. It should be stressed that the leftmost cavity, that is, the "single" 
Fabry-Perot resonator of Figure 4, is initially devoid of any radiation (i.e., it is initially an empty cavity), 
so that the emission of a macroscopic amount of signal microwaves through the left hole from the left 
side of the apparatus, would be a dramatic manifestation of the dynamical Casimir effect, in which the 
observed signal output must have built up exponentially starting solely from vacuum fluctuations inside 
this initially empty resonator. Since the dynamical Casimir effect is closely related to Hawking radiation 
according to lf32ll . an observation of parametric oscillation resulting from the moving SC membrane in 
Figure 4 would be a very interesting result from the point of view of quantum field theory. 



4 The gravitational dynamical Casimir effect, and the generation 
of coherent gravitational radiation 

In this final section, we speculate that the above ideas can be extended to include the case of gravitational 
radiation. The physical concept that ties all these ideas together is the crucial use of the DeWitt minimal 
coupling rule in all of them. 

In particular, we briefly comment on the possibility of extending the "separated parametric oscillator" 
idea for generating EM microwaves by means of the vibrating SC membrane placed inside the extremely 
high Q "triple" SC cavity, as depicted in Figure 4, to the much more speculative idea of generating GR 
microwaves using the same vibrating SC membrane inside the same "triple" SC cavity. This extension 



is based on the fact that the interaction Hamiltonian in (39i is mathematically identical to that 



of the interaction Hamiltonian Ha a m (38 1. Furthermore, we are assuming that it is permissible for 



gravitational radiation fields to be second quantized (see (42 1). 
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However, for this extension of the parametric oscillator idea to work, it is crucial that the walls 
SC cavity, including the surfaces of the moving SC membrane, reflect GR microwaves with as high a 
reflectivity as in the case of EM microwaves. In the paper "Do mirrors for gravitational waves exist?" 
it was predicted that even thin SC films are highly reflective mirrors for GR plane waves. This 



surprising prediction was based on the DeWitt minimal coupling rule (20i applied to the Ginzburg- 
Landau theory of superconductivity. The "off-diagonal long-range order" (ODLRO) l34l nature of the 
Cooper pairs causes these pairs to behave differently from the ions in the ionic lattice, for which ODLRO 
does not exist. As a result, inside the SC thin film, the Cooper pairs, which exhibit constructive AB 
interference, do not undergo geodesic motion, in contrast to the ions, which do undergo geodesic motion, 
in response to incident GR radiation. This difference in the internal motions of the Cooper pairs and of 
the ions inside the SC in the presence of GR radiation, leads to a charge separation effect induced by an 
incoming GR plane wave, such that a huge back-action of the SC film on the GR wave that causes its 
reflection, results. 

If such SC mirrors for GR waves were indeed to exist in Nature, then moving SC mirrors would not 
only be able to do work like a piston on these waves, but would also simultaneously lead to a Doppler 
effect that leads to the exponential amplification of these waves above the threshold for parametric os- 
cillation, as explained above. Thus, a laser-like generation of coherent GR waves starting from vacuum 
fluctuations should become possible. If so, a Hertz -like experiment for GR radiation at microwave fre- 
quencies 1 15 1 would become feasible to perform. 
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